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ABSTRACT 

We investigate chemical evolution in Milky Way-like galaxies based on the cold dark matter model in which 
cosmic structures form via hierarchical merging. We introduce chemical enrichment due to type la supernovae 
(SNe la) into the Mitaka semi-analytic galaxy formation model developed by Nagashima & Yoshii. For the first 
time we derive distributions of stellar metallicities and their ratios in Milky Way-like galaxies treating chemical 
enrichment due to SNe la in a hierarchical galaxy formation model self-consistently. As a first attempt, we 
assume all SNe la to have the same lifetime, and assume instantaneous recycling for type II supernovae (SNe 
II). We find that our model reproduces well the metal abundance ratio [O/Fe] against [Fe/H] and the iron 
metallicity distribution function in the solar neighborhood. This means that the so-called G-dwarf problem is 
resolved by the hierarchical formation of galaxies, and a gas infall term introduced in traditional monolithic 
collapse models to solve this problem is well explained by the mixture of some physical processes such as 
hierarchical merging of dark halos, gas cooling, energy feedback and injection of gas and metals into hot gas 
due to supernovae. Our model predicts more oxygen-enhanced stars in bulges at [Fe/H] ~ than in disks. This 
trend seems to be supported by recent observations while they have still uncertainties. More data in number and 
accuracy will provide independent and important constraints on galaxy formation. For the better understanding 
of the chemical enrichment due to SNe la in hierarchical galaxy formation, we discuss how physical processes 
affect the metal abundance ratio by varying the lifetime of SNe la, the star formation timescale and the strength 
of supernova feedback. We find that the supernova feedback plays a key role among them and therefore there 
is no one-to-one correspondence between star formation histories and stellar metallicity-ratio distributions. 

Subject headings: stars: abundances - Galaxy: abundances - Galaxy: evolution - Galaxy: formation - large- 
scale structure of universe 



1. INTRODUCTION 

The chemical compositions of stars and gas, as is widely 
known, provide important clues to understanding the forma- 
tion and evolution of galaxies. In particular the abundance 
ratio of a-elements such as oxygen to iron is a useful probe of 
the formation history of galaxies. The main sources of met- 
als are considered to be supernovae (SNe), especially type II 
SNe (SNe II) and type la SNe (SNe la). They have different 
abundance patterns, the former produces a-elements as well 
as iron and the latter produces hardly any a-elements. There- 
fore, the measurement of [a/Fe] is closely related to the ratio 
of SNe II to SNe la. While the progenitors of SNe II are mas- 
sive stars larger than ~ 10M©, those of SNe la are considered 
to be white dwarfs in binary systems. Therefore the explosion 
rate of SNe II is almost proportional to the star formation rate 
of their progenitor stars at that time, but there is a time lag 
between SNe la explosions and the formation of their progen- 
itor stars. Since the evolution timescales of those elements 
are independent and different, we can obtain useful informa- 
tion on galaxy formation and star formation histories from the 
abundance ratio of those elements. 

Because SNe II explode almost instantaneously, the 
chemical evolution process has usually been formulated 
by assuming that massive stars immediately explode right 
after their format ion (the instantan eous recycling [IR] 
apDi-oximation) jTinsle\l 19801 lLacev & FarHfT983l fl985l) . In 
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those works, galaxies are usually assumed to grow up in 
monolithic clouds, in which timescales of gas infall and 
star formation (SF) are typically treated as free parameters. 
Since it is possible to observe and resolve individual stars 
in the solar-neighborhood, many studies have focused on 
evolution of the galactic disk. By adjusting those parame- 
ters an d by inco rporatin g stellar p opulation synthesis tech- 
niques ( Arimoto & Yoshi ii 1 986ll 1 9871) . detailed galaxy evolu- 
tion models that reproduce observa tions have been developed 
(Arimoto, Yoshii & Takahara 1992). 

In the framework of such traditional galaxy evolution 
models, several authors explicitly considered chem- 
ical e nrichment due to SNe la JGre ggio & Renzini 
1983 1 : IMatteucci & Greggio| 119861 IMatteucci & Francoisl 
19891 iTsuiimoto et alJ 119951 iPagel & Tautvaisiend 119951 
Yoshii, Tsuiimoto & Nomoto 1996). They included de- 
layed production of metals due to SNe la in the galaxy 
evolution models. By comparing their models with so- 
lar neighborhood stars in the metallicity distribution 
function (MDF) and abundance ratio su ch as [O/Fel 
against metallicity. | Pa flel & Tautvaisiene ( 1995) and 
Yoshii, Tsuiimoto & Nomoto ( 1996) inferre d a typical life- 
time o f SNe la of 1.3~1.5 Gyr. Although Kobavas hi et alJ 
( 1998) has clai med that SNe la s t art to explode at the age of 
about 0.6 Gvr. lGal-Yam & Mao3j2004l) recently derived that 
the lifetime should be longer than about 1 Gyr at the 99 per 
cent confidence level by comparing the cosmic star formation 
history with the cosmic SNe la rate. 

Considered from the viewpoint of cosmological structure 
formation, however, these models are still only phenomeno- 
logical. Recent theoretical and observational studies of struc- 
ture formation have revealed that the universe is dominated 
in mass by cold dark matter (CDM) and that baryonic ob- 
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jects such as galaxies form in virialized objects of dark mat- 
ter called dark halos. In the CDM model, since smaller-scale 
density fluctuations in the early universe have larger fluctua- 
tion amplitudes, larger halos form through mergers of smaller 
dark halos (the hierarchical clustering scenario). This sug- 
gests that large galactic gas clouds would have formed at low 
redshift by mergers of subgalactic clumps which have formed 
at higher redshift. Therefore we may have to modify our pic- 
ture of chemical evolution in galaxies as well as that of galax- 
ies themselves. 

Recently, based on a semi-analytic (SA) approach, galaxy 
formation models in the hierarchical clustering scenario have 
been developed, in which the formation histories of individual 
dark halos are followed by using a Monte Carlo method and 
physical processes such as gas cooling and star formation are 
taken into account in the histories of dark halos. This is a 
natural consequence of adopting the CDM model because we 
are not able to freely assume formation histories of dark halos 
as far as we would like to construct galaxy formation models 
consistent with the CDM model. Many authors have found 
that such SA models reproduce well various characteristics of 
galaxies at the present and at high redshift such as luminosity 
functions, gas fracti ons, size distributi ons, and faint galaxy 
number co unts (e.g..lKauffmann. White & G uiderdonil|1993|: 
Cole et alJU994l 120001 INagashim a. Gouda & Sugiuralp99| 
Nagashima et alJ l200ll 120021 iNagashima & Yoshiil I20041 
Somerville & Primack 1999; Somerville. Primack & Fabei 
2001). In recent SA models, chemical enrichment is con- 
sidered, but generall y only with th e IR approximation. 
iKauffmann & CharloJ (119981) and iNagashi ma & Goudal 
d200ll) investigated color-magnitude and metallicity- 
magnitude relations of cluster elliptical galaxies. 
Thi s is extended to dw a rf sph eroid als, Mb ~ —10 , 
bv INagashima & Yoshiil J2004I) iKauffmannl ill 9961) 
Some rville. Primack & Faberl (120011) and lOkoshi et alJ (12004 
consider chemical evolution in spiral galaxies, a part of which 
should be identified as damped Ly-a systems. Some of them 
found good agreement with observations. 

Pioneering work taking into account chemical en richment 
due to SNe la in a SA model was carr ied out by Thomas 
(fT999h and iThomas & Kauffmannl (f7999h . In those papers, 
they picked out averaged and individual formation histories 
of dark halos, and then, assuming the closed-box chemical 
evolution, they followed star formation and chemical enrich- 
ment histories. Thus, while they took into account merg- 
ing histories of galaxies, the outflow of gas caused by SNe 
or the supernova (SN) feedback that has been realized to be 
an primarily important pr ocess in galaxy formation, wa s not 
considered. As shown bvjKa uffmann & Chariot ( 1998|) and 
Nagashima & Gouda] (120011) . SN feedback significantly af- 
fects chemical enrichment, at least due to SNe II. This sug- 
gests that the closed-box model has limitations in analyzing 
chemical evolution in a realistic situation. 

As a complementary approach to SA modeling, chemo- 
dynamical simulations including SNe la have been de veloped 
(e.g..lRaiteri. Villata & Navarrj 19961 iBerczildl 19991 iKawatal 
1 2001 1 : iLia. Portinari & Carrarol I2002t iNakasato & Nomotol 
120031 lKobavashill2003HOkamoto et alll2005l) . Although this 
approach has an advantage in resolving spatial structure, be- 
cause of the limitation of numerical resolution and numeri- 
cal tech niques themselves, fu rther improvements are still re- 
quired ( Okam oto et alJ 12003). On the other hand, the SA 
model is free from numerical effects and limitations. This is a 
great advantage for understanding global properties of galaxy 



formation because in principle the SA model can investigate 
from dwarf galaxies to galaxy clusters simultaneously. 

In this paper, we construct a fully self-consistent treatment 
of chemical enrichment with a SA mode l. The basic model i s 
the Mitaka model presented by Nagash ima & Yoshiil |2004), 
which includes a Monte Carlo realization of the merging 
histories of dark halos, radiative gas cooling, quiescent star 
formation and starbursts, mergers of galaxies, chemical en- 
richment assuming the IR approximation, size estimation of 
galaxies taking into account the dynamical response to gas re- 
moval during starbursts and stellar population synthesis. As a 
first attempt, we assume that all SNe la have the same lifetime. 
This is, of course, a rather simplified model. In reality the life- 
times of SNe la are considered to have a broad distribution, 
and there is ev en a claim that low-m etallicity environments 
inhibit SNe la (Kobavashi et a Ol998l) . Our model, however, 
has an essential characteristic of SNe la, that is, the time lag 
between star formation and the explosion. It enables us to see 
how such a delayed explosion of SNe la affects the chemi- 
cal enrichment and abundance pattern in galaxies. Using this 
model, we focus on the chemical enrichment in Milky Way 
(MW)-like galaxies. This has a particularly meaning because 
there are many non-trivial effects on the chemical enrichment 
in the hierarchical formation of galaxies. Massive galaxies 
such as the MW form not only via gas cooling but via merg- 
ers of pre-galactic sub-clumps with stronger efficiencies of the 
SN feedback than those in massive galaxies residing in deeper 
gravitational potential wells. Thus, the main purpose of this 
paper is t o compare the SA model with a mo nolithic collapse 
model of Yoshii, Tsuiim oto & Nomotol ( 119961) . This will en- 
able us to see how the hierarchical formation process affects 
galactic metal enrichment due to SNe la based on the CDM 
model. Further comparison to other objects will be done in 
subsequent papers. 

This paper is outlined as follows. In §2 we describe our SA 
model. In §3 we provide a detailed prescription of the chem- 
ical enrichment due to SNe la. In §4 we show the luminosity 
function of galaxies that should be compared with that in the 
Local Group, and properties of MW-like galaxies. In §5 we 
compare our model with observations in the [0/Fe]-[Fe/H] 
plane and the [Fe/H] distribution in a statistical sense. In §6 
we briefly discuss individual galaxies. In §7 we investigate 
parameter dependences of the main results. In §8 we provide 
a summary and conclusions. 

2. SEMI- ANALYTIC MODEL 

Here we briefly describe how to form galaxi es in our m odel, 
which is based on the Mitaka model (Nagashi ma & Yoshiil 
2004). In a CDM universe, dark matter halos cluster gravi- 
tationally and grow in mass through their mergers, depending 
on the adopted power spectrum of the initial density fluctua- 
tions. In each of the merged dark halos, radiative gas cooling, 
star formation, and gas reheating by supernovae occur. The 
cooled dense gas and stars constitute galaxies. These galaxies 
sometimes merge together in a common dark halo, and then 
more massive galaxies form. During these processes, chem- 
ical enrichment due to SNe II is solved by assuming the IR 
approximation, and the enrichment due to SNe la by using 
the past star formation histories of individual galaxies (see 
the next section). Repeating these processes, galaxies form 
and evolve to the present epoch. The Mitaka model repro- 
duces many observations such as luminosity functions, cold 
gas fraction, sizes of disks and spheroidals, surface bright- 
nesses, the Faber-Jackson relation, mass-to-light ratio and 
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faint galaxy number counts. The details , apart from chemical 
enrich ment due to SNe la are found in Nagashi ma & YoshiH 
<200|). 

Throughout this paper, we consider only a A-dominated 
CDM (ACDM) cosmology with fl = 0.3, il A = 0.7, h = 
0.7,(78 = 0.9 and Hh = 0.02/i -2 , and a power spectrum 
given by Sugiv amal (I1995I) . which takes into account the ef- 
fects of baryon. 

The merging his tories of dark halos are calcu lated by using 
a method given by Somerville & Kolatt ( 1999). The timestep 
is Az = 0.06(1 + z), which corresponds to the dynamical 
timescale of dark halos collapsing at z. Only dark halos with 
circular velocity larger than V\ ow — 40 km s _1 are regarded 
as isolated halos. While the value of V\ ow hardly affects the 
chemical evolution in MW-like galaxies that we consider in 
this paper, we will compare the luminosity function of galax- 
ies in the Local Group (LG) with the model luminosity func- 
tions with various Vi ow in Section^ The circular velocities of 
root halos, which are halos at z = 0, is set to be V C y TC = 220 
km s _1 , which is nearly equal to the rotation speed of the 
MW. 

The cycle among baryonic components in dark halos is 
computed as follows. Diffuse baryonic gas is shock-heated 
to the virial temperature estimated from the depth of the grav- 
itational potential well of its host dark halo when it collapses. 
This is called hot gas, whose mass is A/hot- The hot gas dis- 
sipates its energy and cools by radiation only for halos with 
Kirc < V cu t = 250 km s _1 , to avoid formation of unphysi- 
cally giant galaxies, in each timestep. The cooling rate is es- 
timated by using the metallicit y-depe ndent cooling functions 
given by Sutherland & Dopita ( 1993). This cooled gas con- 
centrates in the vicinity of the center of the dark halo and 
constitutes a galactic disk. The mass of cold gas is M^id- 
Disk stars form by consuming the cold gas with a star for- 
mation timescale t* . In proportion to the star formation rate 
ip = M co id/T*, massive stars explode as SNe II. The energy 
released by SNe II reheats the cold gas and a part of the cold 
gas is absorbed into the hot gas (SN feedback). The reheating 
rate is given by M re heat = Pi*, where = (Vdisk/Vhot)~" hot 
and Vhot and ahot are f ree parameters. In this paper the star 
formation timescale is assumed to be t°(1 + 0), where r° 
is a free parameter, which is calle d the constan t star forma- 
tion model (CSF) in Nagashi ma & Yoshiil (120041) . Because (3 
is much larger than unity in low mass halos owing to their 
shallower gravitational potential wells, most of cold gas is re- 
heated and only a small fraction of it turns into stars. The 
above processes are computed during each timestep. 

When a dark halo collapses by merging two or more dark 
halos, the hot gas components in the progenitor halos are im- 
mediately merged and constitute a hot gas component in the 
new dark halo. In contrast, galaxies do not merge together 
immediately. At first we define a central galaxy in the newly 
collapsing dark halo as the central galaxy in the most massive 
progenitor halo, and the rest of the galaxies are regarded as 
satellite galaxies. Two modes of galaxy mergers are consid- 
ered. One is a merger between a central and satellite galaxies, 
in which a satellite falls into a central galaxy after los ing its 
energy by dynamical friction ( Binn ev & Tremainell987l) . An- 
other is that between satellite galaxies by random collisions 
(Maki no & HutlH997l) . We estimate the timescales of the dy- 
namical friction and the random collisions for all satellites. 
Satellites with sufficiently short merger timescales compared 
with the timescale of subsequent collapse of the dark halo 



TABLE 1 
Model parameters 



Vhot a hot 




yn,o 


J/II,Fe 




150 km s" 1 4 


1.3 Gyr 


7.19 x icr 3 


3.40 x 10" 4 


6.67 X 10" 4 



NOTE. — Vhot, a hot and is determined by matching the luminosity functions 
of galax ies and the mass fraction of cold gas in spiral galaxies with observations. See 
Nagashima & Yoshii ■ 2004 ) for details. The above chemical yields are the same as those adopted 
by Yoshii, Tsuiimoto & Nomoto i 1996). In the case of Salpeter's IMF ISalpeteJ 195^1 with a 
slope of 1.35, according to their paper, the lower and upper mass limits correspond to 0. 03 Mq 
and 50Mq, respectively. The so-called "A" parameter determining binary fractions is 0.015- 
0.055. 



merge with a central or other satellite galaxies. If the mass 
ratio of stars and cold gas of a smaller galaxy to a larger one, 
/ = m sma u/mi arge , is larger than /bulge, a starburst occurs 
with SN feedback and all stars go into the bulge component. 
Then no cold gas remains. Otherwise the smaller satellite is 
simply absorbed into the disk component of the central or 
larger satellite without any activities of star formation. We 
adopt /bulge = 0.4. 

In each timestep, the disk size is estimated as a rotationally 
supported disk under the assumption of specific angular mo- 
mentum conservation for cooling hot gas. Initially, hot gas is 
assumed to have the same specific angular momentum as that 
of its host dark halo. The angular momentum of dark halos 
has a log-normal distribution in terms of the so-called dimen- 
sionless spin parameter. If the estimated disk size is larger 
than before, we renew the size. At the same time, we also 
set the disk rotation velocity to be the same as Vci rc of the 
dark halo. When a major merger occurs, the bulge size is esti- 
mated by assuming energy conservation. Then a starburst oc- 
curs. Because the mass of the galaxy after the starburst is dif- 
ferent from that before, we consider the dynamical response 
to the mass loss of the size and velocity dispersion. To do 
this, we use the formalism of the dynamical response taking 
int o account the underly ing dark matter potential developed 
by Nagashima & Yoshii (2003). This process is required to 
reproduce observed characteristics of dwarf spheroidals such 
as surface brightness and velocity dispersion. 

Through the above processes, we obtain star formation his- 
tories of individual galaxies. From those we estimate the 
luminosities of galaxies using a stellar population synthesis 
techni que. We adopt the simple stellar populations given by 
Kodama & Arimotol ill 997ft . 

These p rocedures are the same as those in 
Nagashima & Yoshii (2004). Readers who are inter- 
ested in the details will find them in that paper. The key 
parameters of the model are tabulated in Tabled 

3. CHEMICAL ENRICHMENT 

The chemical enrichment i s treated in a similar manner to 
Yoshii. Tsuiimoto & Nomoto ( 1996), in which they consid- 
ered chemical enrichment of oxygen (O) and iron (Fe) due to 
both SNe II and SNe la using t he infall model in the f rame- 
work of monolithic collapse ( Arimo to. Yoshii & Takahara 
In their paper, the basic equations describing the gas 
fraction, f g , and metallicities for the z-th elements, Zi, are 
written as 

%=-e^(t)+A(t), (1) 
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= -aZi(t)ip(t) + Z A>i A{t) + ay n ,Mt) 



d(Zifg) 

dt 



+ayi a .i / tp{t - ti a )g(ti a )dti a , (2) 
Jo 

where A(t) denotes the gas infall rate, g(ii a ) is a distribution 
function of the lifetime of SNe la, t\ a , a is the locked-up mass 
fraction in low-mass stars and dead stellar remnants and z/n,i 
and yi^i are the chemical yields of the i-th element produced 
by SNe II and SNe la, respectively (see equations 1 and 2 in 
their paper). The star form ation rate ip is assumed to obey the 
Schmidt-law .Schmidtll963l) . 

where v is a rate coefficient. The infall term is assumed to be 

(i/ii„) 7 exp(-Viin) 



A(t) 



(4) 



t in r( 7 + i,T G /i in ) ' 

where t- m is the infall timescale, T(a, b) is the incomplete 
gamma function and Tq is the age of the Galaxy. In 
the followings, we compare our results with their k = 1 
model having k = 1, (av)^ 1 = 3.01 Gyr, t m = 5 
Gyr, 7 = 1, and a constant lifetime of SNe la, g(t\ a ) = 
8 D (ti a - 1.5Gyr). This value of the lifetime of SNe la 
is consis tent with a theoretically estim ated value using a 
model of Lia, Portinari & Carraro (2002), which is based on 
Greggio & Renzini ( 1983), assuming Kennicutt's initial mass 
function (IMF) of stars (Kennicutt 1983). We have estimated 
that the mean lifetime is about 2.4 Gyr and the median about 
1.8 Gyr. 

In our model, because of the SN feedback, these equations 
are modified. During disk star formation or a starburst in a 
timestep that corresponds to Az = 0.06(1 + z), chemical 
enrichment due to SNe II is solved by assuming the IR ap- 
proximation, 



dAL 



cold 



dt 
dM hot 



dt 



d(M t 



cold^cold.-i 



dt 

rf(Afhot-Zhot,i) 



dt 



= -(a + /.)*/;(.), 
= -(a + (3)Z coMA ip(t) 



(5) 



(6) 



otyn,iip, (7) 



(8) 



where Z C oid,i and -Zhot.i are metallicities for i-th element of 
cold and hot gas, respectively. The star formation rate tp is 
given by M co \ A /t^. Equations (|5}-(|8} do not have infall terms 
because we add the mass of hot gas that should cool during 
a timestep Az to the mass of cold gas at the beginning of 
each timestep. The timestep is sufficiently short to describe 
smooth gas accretion. Therefore no infall term is required as 
an approximation. The mixture of merging histories of dark 
halos, gas cooling and SN feedback should correspond to the 
infall process with a long timescale, ~ 5Gyr. The roles of 
the physical processes are clarified in Section 6. In addition 
to this process, metals released from SNe la are put into the 
cold gas. The mass of i-th element released during a timestep 
[tj-utj] is 



AMi = ayi a 



4>(t - t la )dt. 



(9) 



In order to estimate this quantity, star formation histories of 
individual galaxies are stored. Again, we would like to stress 



that the above procedures of chemical enrichment due to both 
SNe la and II are computed in each timestep. 

Since one of the main purposes of this study is to com- 
pare our results with a monolithic cloud model given by 
Yoshii. Tsuiimo to & Nomotol (119961) . we adopt the same 
chemical yields empirically determined by their analysis, 
which are also tabulated in Tabled Thus we also impose 
two nucleosynthesis constraints they adopted: yn,o/yn. Fc = 
const., and yi a ,o = 0. Returned mass from SNe la is sim- 
ply treated together with that from massive stars, that is, it 
is included in a parameter for the instantaneous recycling, 
1 — a. For simplicity we do not take into account the fol- 
lowing effects that might affect chemical evolution of galax- 
ies. One is energy feedba ck due to SNe la, which was 
investigated in detail in lKawata & Gibson! (I2003albl) . Sec- 
ond is inhomogeneous star formation in a gas cloud consid- 
ering propagation of star formation sites, which is another 
way to extend the simple traditional ga l axy evolution model 
llTsuiimoto. Shigevama & Yoshiil 11999. llshimaru & Wanaiol 
1999). Although it will probably provide a broader dis- 
tribution in metal abundances, we concentrate, in this pa- 
per, on how the hierarchical formation process affects the 
chemical enrichment. Third is the metallicity effect that a 
low metallicity enviro nment inhibits the explosion of SNe la 
( Kobavashi et al l 19981) . which might affect the estimation of 
t\ a . These effects would be worth considering in future. 

4. LUMINOSITY FUNCTION AND GLOBAL CHARACTERISTICS OF 
MW-LIKE GALAXIES 

The model we use h ere is the same as that in 
Nagashima & Yoshiij ( 120041) except for the chemical enrich- 
ment process, which do not affect any results in that paper. 
Therefore we do not repeat showing the same figures as shown 
in that paper. 

In the following, we pick out galaxies in dark halos of 
Vdrc = 220 km s _1 whose central galaxies are identified 
as spiral galaxies wit h i?-band bulge-to-total lum inosity ra- 
tio smaller than 0.4 (iBaugh. Cole & Frenkl 1 19961) and with 
-22 1 < Mi - 51ogfe < -21.6 (ISomerville & Primackl 
119991) and disk rotatio n velocities 210 < T^sk/km s _1 < 230 
(Benson et al. 2002a). Hereafter, we call a dark halo satisfy- 
ing the above criteria a LG-halo. We realize 100 LG-halos 
and average over the realization. 

Before entering the details, we see the characteristics of 
MW-like galaxies and their satellites. In Figure^ we show 
a luminosity function of galaxies in LG-halos. The solid line 
indicates the model luminosity function and the histogram the 
observed one compiled by M. Irwin 4 . Because the observed 
data are for the MW and M3 1 groups, luminosity functions for 
models are simply doubled to be compared with the observed 
data. 

Our model shows a broad agreement with the observed lu- 
minosity function in the LG at Mb < —15, although at fainter 
magnitude too many galaxies are predicted. While the ob- 
served luminosity functions actually depend on the definition 
of the Local Group, our conclusion about the chemical evolu- 
tion of the MW-like galaxies does n ot change as shown be low. 
The detailed discussion is given bv lBenson et alJ ( 2002b). 

Such overabundance of low luminosity galaxies is easily 
eliminated by taking into account some physical or observa- 
tional effects as shown by many authors. One of the important 
effects is photoionization due to ultraviolet (UV) background 

4 See http://www.ast.cam.ac.uk/~mike/local_members.html 



SNe la IN A HIERARCHICAL GALAXY FORMATION MODEL 



5 



si 
v 



60 

40 

20 


60 

40 

20 


60 

40 

20 



model 
Obs. (Irwin 




FIG. 1. — Cumulative Luminosity function in the LG. (a) The solid line 
indicates the reference model. The histogram denotes observed luminosity 
function in the LG compiled by M. Irwin, (b) Dependence on Vi ow . The solid 
line is the same as model in (a), the reference model of Vi ow = 40 km s — 1 . 
The dashed, dot-dashed and dotted lines indicate variants of the reference 
model with V| ow = 30, 50 and 70 km s -1 . (c) Dependence on S t h ■ The 
solid line is the same as model in (a), the reference model of 5th = °°- The 
dashed, dot-dashed and dotted lines indicate variants of the reference model 
with 5th = 26, 25 and 24 mag arcsec -2 , respectively. 



or reionization (e.g.. [N agashima. Gouda & Sugiural 119991 
iNagashima & Gouda 2001; Somerville 2002; Benson et al.l 
2002a b, 2003). Here we show the effect of the Jeans mass 
increasing due to reionization in a simple way by changing 
Flow Note that in reality this effect should turn on only af- 
ter the ep och of reionization. Indeed, this simple model is 
similar to Kauffmann. White & Guiderdoni ( 1993), in which 
they assumed that gas cooling is prohibited in dark halos with 
Kirc < 150 km s -1 at redshifts between 1.5 and 5 to be 
consistent with the LG luminosity function in the standard 
CDM model. In Figure ^i, we show luminosity functions 
with different V\ ow . As shown in the panel the solid, dashed, 
dot-dashed and dotted lines indicate Vi ow = 40 (the reference 
model), 30, 50 and 70 km s -1 , respectively. Evidently large 
Viow inhibits the formation of dwarf gal axies. Ano ther effect 
is inco mpleteness in the observations. Nagashi ma & Yoshill 
(2004) predict that there are many dwarf galaxies with sur- 
face brightnesses too low to detect. In the bottom panel, Fig- 
ure^, we show luminosity functions with different threshold 
surface brightnesses, 5th, which determines the lowest sur- 
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FIG. 2. — Characteristics of MW-like galaxies. The crosses indicate each 
MW-like galaxy, (a) Disk rotation velocity Vdisk- (b) Cold gas fraction, (c) 
Effective radius of disks. 



face brightness of galaxies to be detected. The solid, dashed, 
dot-dashed and dotted lines indicate Sth = oo, 26, 25 and 
24 mag arcsec -2 , respectively, where surface brightness is 
defined as that within the effective ra dius. Altho ugh the 
lowest one is about 26 mag arcsec -2 ( Mateol [l998t) . clearly 
observations are not a surface bri ghtness-complete sample. 
As shown in Nagash ima & Yoshill J2004I) . the lowest surface 
brightness of model galaxies reaches about 30 mag arcsec -2 
in B-band at Mb — —10. In addition, some of these low 
surface brightness galaxies might be disrupted by tidal inter- 
action with other galaxies (e.g., Martin et al. 2004). These fig- 
ures suggest that the so-called satellite-galaxy problem might 
not be so serious, at least for luminous objects. Anyway, 
changing V\ ovl and 5th hardly affects our results on chemical 
enrichment of the MW-like galaxies. 

In Figure|2] we show the disk rotation velocity, Vdj s k, cold 
gas fraction, / gas , and effective radius of the disk, rdisk, of the 
MW-like galaxies. The crosses denote individual MW-like 
galaxies. Since the disk rotation velocity is adjusted to be the 
circular velocity of the host dark halo only when the mass of 
disk increases, it is not always the same as y c i rc , which is cho- 
sen to be 220 km s -1 at z = 0. The rotation velocity obtained 
is distributed around the value 220 km s -1 that is similar to 
the observed rotation velocity of the Galactic disk. The cold 
gas fraction is about 10-15% in mass, and the effective radius 
is around 7/i -1 kpc. These are very similar to the observed 
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properties of the MW. Thus our MW-like galaxies are good 
probes enough to investigate the chemical enrichment in the 
MW. 

5. COMPARISON WITH SOLAR NEIGHBORHOOD AND BULGE 
STARS 

In this model we treat bulge and disk stars separately, so 
hereafter we show results for both separately, keeping in mind 
there are still large uncertainties in observation of bulge stars. 
In the followings, we compare our resul ts not only with a the- 
oretica l infall model given by lYoshii. T suiimoto & Nomotol 
(1996) but also with stars in our Galaxy. We focus on the 
difference between the hierarchical and monolithic models, 
rather than comparisons with observations. We assume that 
our Galaxy is a representative, typical sample of galaxies of 
MW-like galaxies under consideration. Future observations 
will make clear whether this assumption is valid or not. 

Firstly, in this section, we show mean properties of MW- 
like galaxies. Properties of individual galaxies are discussed 
in f|6] 

In Figure [3^ we show the distribution of disk stars in the 
[0/Fe]-[Fe/H] plane by contours. By tracking star formation 
histories and metallicity evolutions, we obtain distributions of 
stars in the plane for individual galaxies. The figure shows 
the average of 100 MW-like galaxies. The levels of contours 
drawn by the thin solid and dashed lines indicate 0.5, 0.4, 0.3, 
0.2 and 0.1, and 0.02 and 0.005 times the largest number of 
stars in grids, respec tively. The solid lines indicate th e k = 1 
model given by Yoshii, Tsuiimoto & Nomoto ( 1996) (see Ta- 
bles 1 and 2 in their paper for details). The crosses and filled 
tri angles are observed data for solar-neighborhood star s give n 
bv lMelendez & Barbuvl d2002b and Edv ardsson et alJ i 19931) . 
respectively. As is clearly shown in this panel, our model 
reproduce well the observations of disk stars including their 
dispersion. The apparent discrepancy at [Fe/H]< —2.5 is 
caused by the first nucleosynthesis constraint, yn.o/yu,Fc = 
const, for all SNe II (see Section Qj. At the break point the 
model prediction gives slightly higher [O/Fe] than observa- 
tions. Probably if the assumption of a constant t\ & were re- 
laxed, the abundanc e ratio at there would de crease as sug- 
gested by Yoshii. Tsuiimoto & Nomoto ( 1996). The disper- 
sion for the model comes from dispersion in a single galaxy 
and that for different samples. In the next section we show the 
former effect by seeing some individual galaxies. 

It should be noted that a small radial abundance gradients 
of metals in the galactic disk are suggested by observations, 
~ —0.0 85 dex kpc" 1 for iron and ~ —0.07 dex kpc -1 for 
oxygen (Maciel 2001). This might cause a systematic differ- 
ence in the metallicity distribution when averaged over the 
whole disk. In that case, we will need to change the val- 
ues of chemical yields, which are e mpirically determined by 
Yoshii. Tsuiimoto & Nomoto ( 1996) to produce an agreement 
with the observations for solar-neighborhood stars. 

Bulge stars are also distributed like disk stars in the [O/Fe]- 
[Fe/H] plane but the range is slightly different as shown in 
Figure |3ji. In particular, although the fraction is low, there 
are many oxygen-enhanced star s at [Fe/H]~ 0. Interest- 
ingly, [Rich & Mc William (120001) found that [O/Fe] for bulge 
stars, shown by filled inverse triangles, has a broadly simi- 
lar trend with [Fe/H] to disk stars, but some of them show 
oxygen-enhancement, [0/Fe]~ 0.3 even at [Fe/H]~ 0. Re- 
ce nt high-resolution spectroscopic observatio ns of bulge stars 
by Mc William, Rich & Smecker-Hane (2003) reveal a simi- 
lar trend of [Mn/Fe] to that for disk stars. While they ob- 



served less than 10 stars, their observations may support 
our model. Clearly more samples will provide important 

constrai nts on bulge formatio n. Using monolithic collapse 

models, Matteucci & Brocato ( 1990) predicted more oxygen- 
enhanced bulge stars. To obtain this, they assumed higher star 
formation rates and shorter infall timescales for bulge stars 
than those for disk stars. In our model, some bulge stars are 
also formed with a very short timescale, that is, starbursts. 
The difference we predict between disk and bulge stars, how- 
ever, is smaller presumably because not all bulge stars are 
formed by starbursts and some have been formed in progen- 
itor galaxies. Another possibility to form oxygen-enhanced 
stars is a variation of the IMF. Recently the IMF for starbursts 
is suggested to b e top-heavy compared to the IMF for t he disk 
star formation (Bau grTet al.l2005t lNagashim a et al 12005 albl) . 
which means that much more SNe II must explode than ex- 
pected for Salpeter-like IMFs. Because a-elements from SNe 
II are likely to be produced more, metal-rich bulge stars will 
tend to be placed at [0/Fe]> 0. 

In Figures |3j; and d, we show the iron MDF for disk and 
bulge stars, respectively. The solid lines denote the model pre- 
diction. The filled squares and crosses in the panel (c) and the 
filled squares, crosses and pluses in the pane l (d) indicate ob- 
ser ved [Fe/H] distribu t ions g iven bv lWvs e & Gilmore ( 1995) 
and lEdva rdsson e t al.l (1199 3) for di sk stars and IZoc cali et all 
(2003), Ramire z et alJ ll2000l) and ISadler. Rich & Terndrupl 
( 1996) for bulge stars, respectively. All of these distribu- 
tions are normalized for the integrate over the whole range 
to be unity. The distribution for disk stars shows an excel- 
lent agreement with the observed distributions. Agreements 
in abundance ratio and MDFs are quite encouraging to model 
the formation of the galactic disk in the framework of SA 
models. For bulge stars, there remains discrepancy among 
observations the mselves. O ur model i s consistent only with 
the data from Zoccali et al. (2003) and Ramirez et al. (2000). 
The data from Sadler. Rich & Terndrupl ill 9961) shows more 
metal-rich stars than the other two and our model. 

While the predicted MDF for bulge stars is rather similar to 
that for disk stars, the width is clearly wider than that of disk 
stars. The observations of bulge stars show a similar tendency 
to the model, although there seem to be still uncertainties in 
observations outside the range —1 <[Fe/H]< 0. Keeping in 
mind differences between observational data, the model pre- 
diction for bulge stars also agrees well with the observed data. 
In another approach, a recent chemo-dynam i cal sim ulation of 
the galactic bulge by Nakasato & Nomoto (2003) also pre- 
dicts distribution of stars in [O/Fe] -[Fe/H] plane and in num- 
ber against [Fe/H]. In their results, ther e are more oxygen- 
enhanc ed stars than in the observation bv lRich & M c William 
(2000) and our prediction. At the same time, their derived 
[Fe/H] distribution is similar to ours, although the number of 
high [Fe/H] stars is a little larger. Thus, although more data 
will be required to constrain the model of bulge formation, 
we can say that both our SA model and the recent numerical 
simulation broadly agree with the current observations. 

Figure |4] shows averaged star formation histories of MW- 
like gala xies indicated by the solid line. T he infall model 
given by Yoshii. Tsuiimoto & Nomoto ( 1996) is indicated by 
the dashed line. Note that the SFR for the infall model is plot- 
ted beyond the cosmic age given by the cosmological param- 
eters we adopt, ~ 13 Gyr, because their model is normalized 
at 15Gyr. The star formation rates are normalized by stellar 
masses at 13Gyr. SFRs plotted are those averaged over each 
timestep (AM, /a At) and over all realization. Clearly galax- 
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FIG. 3. — (a) [O/Fe] distribution against [Fe/H] for disk stars. The levels of contours drawn by the thin solid and dashed lines indicate 0.5, 0.4, 0.3, 0.2 and 0.1, 
and 0.02 and 0.005 times the largest numbe r of stars in grids , respectively. The solid curve indicates a chemical enrichment model based on a monolithic cloud 
coll apse model give n by Yoshii, Tsuiimoto & Nomoto 1 1996). The circles with error bars and filled triangles denote observations by Melendez & Barbuv 1 2002 ) 
and Edvardsson et al. 1 1993), respectively, (b) [O/Fe] distribution against [Fe/H] for bulge stars. The levels of contours and the solid curve are the same as the 
panel (a). The filled inverse triangles denote observational data by Rich & McWilliam 1 2003). (c) Ir on MDF for disk stars. The s olid line indicates the model 
prediction for disk stars. The filled squares and crosses denote observations by Wvse & Gilmore 1 1995) and Edvardsson et al. 1 1993), respectively, (d) Iron MDF 
for bulge stars. The solid line indicates the model prediction for bulge stars. The filled squares, crosses and pluses denote observations by Zoccali et al. i2003), 
Ramirez et al. |2000) and Sadler, Rich & Terndrut) 1 1996), respectively. 



ies in our SA model have similar star formation histories to 
the monolithic collapse model's. Hence the slow star forma- 
tion in the infall model, which is brought by long infall and 
star formation timescales chosen by obtaining an agreement 
with observations, is naturally explained by the mixture of 
physical processes including merging histories of dark halos, 
gas cooling and SN feedback. While a history of each MW- 
like galaxy has a variety around the averaged history, as we 
will show in the next section, we have found that predicted 
metallicity distributions of disk stars are quite similar to the 
averaged one shown in Figure|3] Since bulge components are 
formed by major mergers that are rare events, it should be 
noted that the results on the bulge stars as just an average. Fu- 
ture observations of bulge stars in other galaxies will provide 
stronger constraint on galaxy formation models based on the 



CDM cosmology. 

Our results suggest that the formation of spiral galax- 
ies is well approximated by the traditional infall model be- 
cause our model provides similar results to those given by 
Yoshii, Tsuiimoto & Nomotol d!996l). This i s consistent with 
the result by Baugh, Cole & Frenk ( 1996), in which they 
traced formation histories of individual galaxies and found 
similar star formation histories of spiral galaxies to the tra- 
ditional infall model (see upper panels of Figure 3 in their 
paper). Note that recent infall models have become more 
complicated, including multi-accretion of gas, in order to re- 
produce observational results. Moreover, to reproduce radial 
gradients of metallicities, the mod els require SFRs dependent 
on th e radius of disks (e.^.. | Chiappini. Matteucci & Grattonl 
U997llFenner & Gibsonll2003l) . Physically, the radial gradi- 
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FIG. 4. — Star formation histories. The solid line indicates the averaged 
star formation history of MW-like galaxies. The dashed line denotes a history 
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disk stars of individual MW-like galaxies against their last major merger 
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ents should be solved by computing complex dynamics of 
disk formation. This requires numerical simulations of galaxy 
formation following dynamical evolution including the mix- 
ing process of gas and metals. 

Finally we shortly discuss relations between averaged 
metallicities of individual galaxies and their merging histo- 
ries. Figure|5]shows mass-weighted mean metallicities of disk 
stars for O (crosses) and Fe (pluses) against their final major 
merger epochs, z mrg . Correlations of metallicities with z mlg 
are very weak for O and negligible for Fe. Particularly, the Fe- 
metallicity is almost independent of z mrg . The O-metallicity 
only for galaxies which experienced a major merger very re- 
cently is slightly lower or has a larger scatter than others. To 
make the reason clear, however, we need more samples of 
high-z MW-like galaxies. The current number of high-z sam- 
ples is still few for quantitative statistical discussion. At this 
stage, what we can say is that there is no clear correlation 
between the metallicity and the final major merger epoch. 

Figure |6] shows the same as the previous figure but for fi- 
nal major merger epochs for their host halos, Zhaio mrg- Simi- 
larly, we do not see correlations between the O-metallicity and 
Zhalo mrg- The negligible or weaker correlation with Zhaio mrg 
than with z mrg can be reasonable because the galaxy merger 
does not directly reflect the halo merger and is lead by other 
physical mechanisms, namely, dynamical friction and random 
collisions. 

Thus, we conclude that the metallicity depends only weakly 
on the formation history of galaxies and their host halos. In- 
stead, we can see much stronger dependence of the metallicity 
on parameters of, for example, SN feedback later. 

6. INDIVIDUAL GALAXIES 

In the previous section, we showed mean properties of MW- 
like galaxies. We here show results for two samples in order 
to see properties of individual galaxies. 

FigureQshows luminosity functions of galaxies in two LG- 
halos hosting MW-like galaxies as central galaxies of ID- 143 
with z mrg = 1.13, Zhaiomrg = 0.79, (Z /Z .q) = 0.89, 
and (Z Fe /Z Fe Q ) = 0.40, and ID-200 with z mrg = 0.59, 
Zhaiomrg = 0*.50, (Z /Z ,@) = 0.99, and (Z Fc /Z FCi& ) = 
0.47, respectively. The luminosity function in ID- 143 halo 
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FIG. 6. — Mass-weighted metallicities of O (crosses) and Fe (pluses) for 
disk stars of individual galaxies against last major merger epochs of their host 
halos, z halo mrg . 



agrees well with the observed one at My < — 15, and that in 
ID-200 halo broadly agrees but slightly overproduces galax- 
ies at My ~ —17. This slight discrepancy between the mod- 
els and the observations should substantially decrease when 
taking into account some effects such as the selection effects 
shown in We also find that difference between samples is 
not significant. 

Figures [8] and [9] show the same figures as Figure [3] that is, 
the distribution of disk and bulge stars in [0/Fe]-[Fe/H] plane 
and MDFs for [Fe/H], but for individual galaxies of ID-143 
and ID-200, respectively. In order to enhance the contribu- 
tion of stars with small fractions, the contour levels are not 
linear but logarithmic. As easily can be seen, some stars are 
located in a low [O/Fe] region at low [Fe/H], particularly for 
ID-143. Such stars can be brought by smaller sub-clumps 
via minor mergers. Because low-mass galaxies with shallow 
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FIG. 7.— Cumulative luminosity functions in LGs of ID- 143 (solie line) FIG. 9.— The same as Figure|S]but for a MW-like galaxy of ID-200. 

and ID-200 (dashed line). No selection effect is taken into account. The 
histogram denotes the same observed luminosity function shown in FigurefTI 
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FIG. 8. — The same as Figure^]but for a MW-like galaxy of ID- 143. Note 
that, in order to enhance contribution of stars with small fractions, the contour 
levels are equally divided in logarithmically. 



that of ID-200's. Thus it should have been hard for ID-200 
galaxy to accrete small sub-clumps as minor mergers until 
present. 

Similar trails can be seen for bulge stars both of ID- 143 and 
ID-200 galaxies, in contrast to disk stars. Because many stars 
in bulges formed in disks firstly and then they were transferred 
to bulges by major mergers, the trails in bulges should reflect 
minor-merger histories onto disks before formation of bulges, 
or, last major mergers. Therefore these galaxies could have 
enough time to accrete sub-clumps until last major mergers. 

The widths of MDFs for both MW-like galaxies are slightly 
narrower than the averaged one shown in Figure [5] In our 
model we assume that metals are well mixed in a galaxy. In 
reality, however, the mixing of metals in a galaxy would not 
be completed and therefore the MDF would be broadened as 
observed. If we take this effect of inhomogeneous metal en- 
richment into account, the MDF will become wider. Note that 
recent o bservations of more than 10,000 solar-neighborhood 
stars by Nordst r5m et alJ fc004l) have suggested a narrower 
MDF than past observations. Although we should be careful 
in sample selection, it may suggest that the metal mixing is 
efficient for disk stars. This should be addressed in future. 



gravitational potential wells suffer strong SN feedback, the 
metal enrichment does not proceed significantly. Therefore 
SNe la start to explode and lower [O/Fe] when [Fe/H] is still 
low. This effect shall be shown in Figure^]for a stronger SN 
feedback case. Although the number of such stars is small, it 
might be discovered as relics of past minor merger histories. 
Note that, of course, this depends on the details of the process 
of minor mergers. 

In contrast to ID- 143, the galaxy ID-200 does not have such 
low-[0/Fe], low-[Fe/H] disk stars. This is because minor 
mergers did not happen after the last major merger and the 
disk grew in a similar way to the infall model of monolithic 
collapse. Merging histories of these galaxies should provide 
insights for understanding the difference in the stellar distri- 
bution in [O/Fe] -[Fe/H] plane. For example, the last merger 
epoch for this galaxy is z mlg = 0.59, and for ID-143 galaxy 
with many minor mergers z mrg = 1.13, which is earlier than 



7. PARAMETER DEPENDENCE 

Here we investigate how the stellar metallicity depends on 
the fundamental model parameters such as the SN feedback 
and star formation timescale, as well as the lifetime of SNe 
la, in order to clarify how physical processes such as star for- 
mation affect chemical enrichment in galaxies. Note that 
adopted parameters below break the agreement with observa- 
tions. We vary the values of the parameters just to see how 
each physical process affects metal enrichment. 

7.1. Lifetime of SNe la 

Firstly, we vary the lifetime of SNe la, t\ a . In Figures 
1101 and fTTl we show the same model as in Figure [5] but for 
ti a = 0.5 and 3 Gyr, respectively. In the case of the shorter 
ii a , because enrichment due to SNe la begins earlier, the break 
point in the [0/Fe]-[Fe/H] plane moves toward lower [Fe/H]. 
On the other hand, in the longer t\ a case, the iron abundance 
initially increases only due to SNe II, so the break point moves 
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FIG. 10. — The same as Figure |3lbut for ti a = 0.5 Gyr. 
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FIG. 12. — The same as Figurel3lbut for = 2.6 Gyr. 



£ 0.5 




(a) disc stars _ 

o Melendez & Barbuy (2002). 
* Edvardsson et al. (1993) 



(b) bulge stars 

T Rich & Mcffilliam (2000) 




I I I I I I I I I I I I 

(c) disc stars 



1 1 1 1 1 1 1 1 1 1 1 

(d) bulge stars 



Wyse & Gilmore (1995) 
Edvardsson et al (1993) 



Zoccali et al. (2003) 
Ramirez et al. (2000) 
Sadler et al. (1996) 





-3 -2-10 1 -3 -2 -1 

[Fe/H] [Fe/H] 

FIG. 1 1 . — The same as Figure|3]but for ti a = 3 Gyr. 



toward higher [Fe/H]. This dependence is basically the same 
as that shown in analyses using the traditional model. This 
also suggests the similarity of formation of spiral galaxies be- 
tween the hierarchical and traditional models. 

In the shorter ir a case, an interesting feature in abundance 
ratio emerges there. Bulge stars has a kind of U-shaped dis- 
tribution with a minimum at [Fe/H]^ —0.5 in [0/Fe]-[Fe/H] 
plane. In this case, iron enrichment due to SNe la is almost 
completed at very high redshift, at which time a large amount 
of cold gas remains in galaxies. Hence most of the stars exist- 
ing at present form in subsequent starburst due to major merg- 
ers between such gas-rich systems. Because starbursts tends 
to make stars at its late phase with the abundance pattern of 
SNe II, if the system is gas-rich, many oxygen-enhanced iron- 
rich stars form during the starburst. Thus this provides a con- 
straint on the lifetime of SNe la, that is, tr a should be rather 
long, > 1 Gyr. 

The dependence of the [Fe/H] distribution on ir a is not as 
simple as the abundance ratio. These figures show that the 



lifetime tu seems to determine the width of [Fe/H] distribu- 
tion, that is, a shorter ti a makes the [Fe/H] distribution nar- 
rower. It could be related to the epoch of star formation. In 
the shorter t\ a case, many intermediate stars explode as SNe 
la and begin to release metals into the interstellar medium be- 
fore the major epoch of star formation compared with the ref- 
erence model, therefore most of the stars have similar iron 
abundances, which makes the [Fe/H] distribution narrow. In 
contrast, in the longer t\ a case, some of the stars form in iron- 
poor environments before most of the iron is released, and 
some of the stars form in iron-rich environments because at 
late epochs the gas fraction is low and iron is released into 
gas-poor systems. Thus the [Fe/H] distribution has a wider 
width. 

7.2. Star formation timescale 

Next, we investigate the dependence on star formation 
timescale. In Figures [I2land[13l we show models which 
are the same as the reference model, except for having twice 
and half the star formation timescale of the reference model, 
r° = 2.6 and 0.65 Gyr, resp ectively. 

As seen in Figure 112b . in the longer star formation 
timescale case, the break point in [0/Fe]-[Fe/H] plane for 
bulge stars moves towards lower [Fe/H], because the star for- 
mation epoch shifts to lower redshift and then the enrichment 
due to SNe la effectively begins earlier. In the case of a shorter 
star formation timescale, by contrast, the break point moves 
towards higher [Fe/H]. Thus we can say that the position of 
the break point is generally determined by a combination of 
the lifetime of SNe la and the star formation timescale, such 
as iia/r°. In the longer case, there is a broader tail show- 
ing oxygen-enhancement at [Fe/H] ~ than in the reference 
model. This is for the same reason as the case of shorter t\ & . 
The difference from the case with shorter tj a is that the dis- 
tribution is not U-shaped but an inverse triangle, because the 
SNe la continue enrichment at lower redshifts. Disk stars also 
have the same dependence on the star formation timescale, 
but we can clearly see those tendencies for bulge stars. 

The [Fe/H] distribution shows a different tendency from 
the dependence on t\ a . In the case of longer and shorter 
t ( J, the distributions typically but slightly have lower and 
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FIG. 13. — The same as Figure l3lbut for = 0.65 Gyr. 
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FIG. 14. — The same as Figure|3]but for the DSF model. 



higher [Fe/H], respectively, at least for bulge stars. Gener- 
ally a longer star formation timescale results in a larger cold 
gas fraction. Then, with the longer star formation timescale, 
the metallicity becomes lower than in a model with a shorter 
star formation timescale. In addition, such a longer star for- 
mation timescale decreases the amount of stars, which means 
that fewer metals are released from SNe la. Thus the star for- 
mation timescale has a different effect from the lifetime of 
SNe la. 

In the reference model, the star formation timescale is as- 
sumed to be constant with redshift. It is interesting to see 
ho w the red shift depende nce af fects the chemical enrichment. 
In Nagashima & Yoshii (2004), another formulation for the 
star formation timescale is considered; it is proportional to 
the dynamical timescale of galactic disk, which is called the 
dynamical star formation (DSF) timescale model in that pa- 
per. Because we assume the disk rotation velocity to be al- 
most the same as the circular velocity of dark halos, which be- 
comes larger toward higher redshift, the dynamical timescale 



T,°=2.6Gyr 
T,°=0.65Gyr 
DSF(1.7Gyr) 
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FIG. 15. — Star formation histories. The solid, dot-dashed and dotted lines 
indicate models with t® = 2.6Gyr and 0.6 5G yr, a nd a model with DSF, 
respectively, corresponding to Figures ll2ll 1 3l and l 1 41 The dashed line denote 
the infall model. 



becomes shorter toward higher redshift. Figure^Jshows the 
results of adopting the DSF model, in which the star forma- 
tion timescale at present is r° = 1.7 Gyr, so as to produce an 
agreement of cold gas fraction in spiral galaxies with obser- 
vation. We find that the position of the break point in [O/Fe]- 
[Fe/H] plane slightly moves toward higher [Fe/H], and the 
slope of [O/Fe] against [Fe/H] seems to be slightly steeper. 
In the DSF model, since most of stars form at higher redshift 
than in the reference model due to the shorter star formation 
timescale at high redshift, enrichment of iron happens quickly. 
Thus the slope is steeper and the typical [Fe/H] is larger than 
in the reference model, but this effect is very small. 

In Figure ^1 we show star formation histories of those 
three models. The solid, dot-dashed and dotted lines indicate 
models with r°=2.6Gyr and 0.65Gyr, and the DSF model, re- 
sp ectively. The dashed line denotes th e infall model given 
bv lYoshii~T suiimoto & Nomoto ( 1996). The major star for- 
mation epochs of the models of t° = 2.6Gyr and 0.65Gyr 
are later and earlier than the reference model shown in Fig- 
ure |4] Because a star formation timescale at high redshift 
in the DSF model is much shorter than that in the reference 
model, the major star formation epoch moves toward even ear- 
lier epoch than the model of r° = 0.65Gyr. From this figure, 
thus we can reach the same conclusion as analyses by the in- 
fall model, that is, when major star formation finishes before 
iron enrichment due to SNe la becomes efficient, the break 
point in [0/Fe]-[Fe/H] moves toward higher [Fe/H] because 
even SNe II alone sufficiently enrich stars. 

As long as we stay varying only the lifetime of SNe la, ti a , 
and the star formation timescale, r*, our SA model seems 
to provide quite similar conclusions to the classical infall 
model. This apparent coincidence, however, breaks when the 
SN feedback is considered as shown in the next subsection. 

7.3. Supernova feedback 

It has been widely realized that SN feedback significantly 
affects galaxy properties. Strong energy feedback suppresses 
chemical enrichment due to SNe II because star formation is 
suppressed by expelling cold gas, and also most of the met- 
als are expelled into the hot gas (Kauffmann & Chariot 1998; 
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FIG. 18. — Star formation histories. The solid and dot-dashed lines indicate 
models with Vhot = 280 and 70 km s _1 , respectively. The dashed line 
denotes the infall model. 
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FIG. 17. — The same as Figure|5]but for Vhot 



Nagashima & Goud al2001l) . Thus it is worthwhile to investi- 
gate how SN feedback affects the enrichment due to both SNe 
II and SNe la. 

In Figures[n[]andE] we show the same figures as Figure|3] 

but for Vhot = 280 (strong feedback) and 70 (weak feedback) 
km s -1 , respectively, where Vhot determines the typical scale 
at which SN feedback becomes effective, that is, SN feedback 
strongly affects galaxies with Vdisk < Vhot- As clearly shown 
in the bottom panels, the [Fe/H] distributions for both disk 
and bulge stars are strongly affected by SN feedback. Strong 
feedback greatly suppresses formation of metal-rich stars and 
thus the chemical enrichment. In contrast, there are too few 
metal-poor stars in the case of weak SN feedback compared 
with the observed distributions, while the difference from the 
reference model is smaller than that of the strong SN feed- 
back case. It should be noted that our MW-like galaxies have 
Vdisk — 220 km s _1 at low redshift. The values of Vh Q t for 
both the reference model and the weak SN feedback case are 
less than the rotation velocity, while in the strong SN feed- 



back case Vhot is larger than the rotation velocity. That means 
that SN feedback is very efficient even at low redshift in the 
strong feedback case. Therefore the difference of the strong 
SN feedback case from the reference model is larger. 

The stronger the SN feedback is, the lower the [Fe/H] at 
the break point is. Consequently, strong and weak SN feed- 
back seem to have similar effects in the [0/Fe]-[Fe/H] plane 
to those of short and long tj a , and long and short r*. Since 
strong SN feedback suppresses the chemical enrichment due 
to SNe II, the enrichment due to SNe la begins at a lower 
[Fe/H] . These apparently correspond to longer and shorter star 
formation timescales, respectively. 

Figure [^] shows star formation histories of the strong 
(Vi 10 t = 280km s" 1 ) and weak (VW = 70km s" 1 ) feed- 
back models indicated by the solid and dot-dashed lines, re- 
spectively. The dashed line denotes the infall model. The 
strong feedback model shows a similar star formation history 
to that of the long star formation timescale model (the solid 
line in Figure [l5t . Both models have break points in [O/Fe]- 
[Fe/H] plane at slightly lower [Fe/H] than that in the reference 
model. The peak of the [Fe/H] distribution of the strong feed- 
back model, however, moves toward lower [Fe/H] compared 
with the long star formation timescale model. 

The weak feedback model shows a similar star formation 
history to that of the DSF model (the dotted line in Figure 
\151 . The breakpoints of both model move toward high [Fe/H] 
compared with the reference model. Again, the shift to higher 
[Fe/H] in the [Fe/H] distribution of the weak feedback model 
is slightly enhanced compared with the DSF model. 

This is because SN feedback strongly affects chemical en- 
richment. From equations 0-©, we obtain a mass-weighted 
mean stellar metallicity, 



(Z*(t)) = Z° c + 



ay 1 — e x — xe 



a + /3 



1 - e" 



(10) 



where x — (a + /3)t/r^ and Z° is an initial metallicity of gas. 
At the limit i/r* — » oo, the mean metallicity becomes 



<Z,(oo))-i£ + 



ay 



(11) 

a + p 

This shows that a strong SN feedback ((3 ^> 1) suppresses 
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chemical enrichment due to SNe II. The effective chemi- 
cal yield taking into account SN feedback can be written as 
VeS = + @l a )- The effect of SN feedback on chem- 
ical enrichment is the reason why, for example, the strong 
feedback model and the long star formation timescale model 
predict different [Fe/H] distributions in spite of their similar 
star formation histories. Therefore, similar star formation his- 
tories do not always provide similar chemical enrichment his- 
tories. Gas outflow induced by SNe II plays a critical role in 
galaxy evolution. 

8. SUMMARY 

We have explored chemical enrichment due to both SNe II 
and SNe la in Milky Way-like galaxies in the semi-analytic 
galaxy formation model. Our treatment of SNe la is fully 
consistent with the galaxy formation model, that is, we solve 
for the recycling of each element among stars, cold gas and 
hot gas based on a ACDM model. It is important to follow 
the chemical enrichment in the framework of the hierarchical 
galaxy formation because there are many non-trivial effects 
on the chemical enrichment caused not only by mergers of 
galaxies but by varying efficiencies of the SN feedback de- 
pendent on the depth of gravitational potential wells of sub- 
galactic clumps. As a first attempt at constructing such a con- 
sistent model, we have assumed that all SNe la have the same 
lifetime, t\ a , that abundance patterns of metals from SNe la 
and SNe II are always the same, and that massive stars instan- 
taneously explode as SNe II an d release metals. This is a nat- 
ural ex tension of the work by Yoshii, Tsuiimoto & Nomoto 
( 1996), in which they self-consistently modeled the chemical 
enrichment due to SNe la and SNe II. 

We have picked out galaxies in dark halos with V c i rc = 220 
km s _1 and having a similar luminosity to the MW. We have 
found that when we impose fi a =1.5 Gyr, the predictions of 
our model for such MW-like galaxies agree well with obser- 
vations of the chemical composition of solar-neighborhood 
stars both in the stellar distribution in the [0/Fe]-[Fe/H] plane 
and in the iron MDF. We would like to stress that the 
other model parameters such as the star formation timescale 
and SN feedback are the sa me as in the fiducial model of 
Nagashima & Yoshii d2004ll . in which they found that the 
model reproduce well many aspects of observed galaxies, 
such as luminosity functions, cold gas fractions and sizes of 
galaxies for local galaxies, and the surface brightnesses, ve- 
locity dispersions, mass-to-light ratios and metallicities of lo- 
cal dwarf spheroidals. This work also shows that the clas- 
sical G-dwarf problem ([van den Bergh 1962; Schmidtl ll963t 
iPagel & Patchet3 ll975l) is fully resolved in the framework of 
hierarchical formation of galaxies, in which the infall term in- 
troduced in the infall models to avoid the G-dwarf problem is 
naturally explained by the mixture of such physical processes 
as clustering of dark halos, gas cooling and SN feedback. Our 
model passes the new tests, that is, the iron MDF and the 
abundance pattern of metals, taking into account chemical en- 
richment due to SNe la. As shown in Sectional the abundance 
pattern provides an independent constraint on galaxy forma- 
tion. Therefore our results would support the scenario of hi- 
erarchical formation of galaxies. While observational data 
used in this work have been limited to solar-neighborhood and 
bulge stars in the Milky Way, future observations of stars in 



other galaxies will provide statistically better constraints on 
galaxy formation. 

In our model, disk and bulge stars are treated separately. 
While the observations of bulge stars still have uncertainties, 
our model also reproduces well both the distribution of stars in 
the [0/Fe]-[Fe/H] plane and the [Fe/H] distribution. In partic- 
ular, although there are only a few bulge stars whose oxygen 
abundance is observed, our model predicts oxygen-enhanced 
stars at [Fe/H]^ 0, but they are not a dominant fraction. In- 
creasing observational data will provide a strong constraint on 
bulge formation. 

As a by-product, we have obtained the luminosity func- 
tion of galaxies in Local Group-halos. Recent high-resolution 
A-body simulations have predicted many more dwarf scale 
dark halos than observed satellite gala xies, and this has 
been consid ered to be a serious problem (Klvpin et al.ll 19991 
Moore et a l .119991) . Such an overabundance problem has been 
solved in the framework of the SA models taking into ac- 
count effects of reionization and incompleteness due to signif- 
icantly low surface brightness ( Somervill el2002tlBenson et"afl 
2002b). Our model also shows an ability to solve this problem 
in a similar manner to those works. 

Th e chemical yields we used ar e the same as those 
in Yoshii, Tsuiimoto & Nomoto ( 1996), in which 
an infall model of mon olithic collapse was used 
( Arimoto. Yoshii & Takahara 1992). The similarity be- 
tween our hierarchical model and the monolithic cloud 
collapse model suggests that spiral galaxies in a hierarchical 
universe should have a similar formation history to that 
m odeled by the monolithic clo ud collapse model, as shown 
in Baugh, Cole & Frenk ( 1996). Star formation histories we 
computed also suggests the similarities. It does, however, 
not always mean that similar star formation histories provide 
similar chemical enrichment histories. As shown in Section 
6.3 and more directly by equation (lilt . SN feedback plays a 
critical role in chemical enrichment. 

In this paper we have concentrated on investigating the 
[O/Fe] relation and the [Fe/H] distribution. Since recent anal- 
yses of solar-neighborhood stars have revealed that the age- 
metallicity relation has a large scatter, it will be useful as a 
next step to see whether our model can reproduce this rela - 
tion as well as the scatter in it (Ibukivama & Arimoto 2002). 
Furthermore, while the analysis in this paper focused only 
on MW-like galaxies, our SA model has the potential abil- 
ity to investigate other systems simultaneously. To investigate 
abundance ratios in such systems as the intracluster medium, 
the stars composing elliptical galaxies, and damped Ly-a sys- 
tems will provide independent, important clues to understand- 
ing galaxy formation. The cosmic explosion rate of SNe la 
will also give a new insight into both galaxy formation and 
observational cosmology. We will pursue these topics in fu- 
ture papers. 
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